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In this paper, the time-dependent coupled-power equations originally 
derived by Marcuse from intuitive arguments are rederived two ways, one 
using the coupled-line equations with perturbation theory and the other 
using the Kronecker product approach of Rowe and Young.* 

I. INTRODUCTION 

Suppose a multimode fiber guide with random mode coupling is 
excited by an optical source at one end at time t = 0. 

This input excitation will produce, at time t, a response of optical 
power in each mode v at a position z from the input end given by 
s v (t, z). 

Marcuse 2 has suggested from intuitive arguments that s v (t, z) should 
satisfy the following differential equation 

£ <«.(*, z)) + -r| <«.(*, i)) - L \K„\*FMt, z)) 

- (L \K^F)Mt,z)), (1) 

where K^ v is the coupling coefficient in the coupled-line equations 
(described below) between modes n and v, F is the spectral height of 
the mechanical perturbation in the fiber geometry responsible for the 
coupling, (s„(£, z)) is the average* power in mode \x at position z, and 
C is the group velocity in mode v. 

In this paper, we derive eq. (1) from the coupled-line equations first 
using perturbation theory and then using the Kronecker product 
approach of Rowe and Young. 3 

The importance of this work is to show that the intuitive eq. (1) 
does in fact follow directly from (i) the coupled-line equations, (ii) 



* Subsequent to the writing of this paper, it became known to the author that many 
of the results included in the paper were independently and concurrently derived by 
R. Steinberg of Columbia University. 1 

* Average for an ensemble of guides with similar gross properties and similar 
excitation. 
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voltage linearity of the guide (from Maxwell's equations), and (Hi) 
ensemble averaging over the random coupling perturbation. It is 
encouraging to note that skillfully framed power-flow arguments such 
as those used by Marcuse lead to the same results as the more cumber- 
some approaches that start from the coupled-line equations. 

It should be pointed out again that, however it is derived, eq. (1) 
describes the average of the flow of power for an ensemble of guides. 
How the flow of power in a particular guide compares to the average 
flow is still an open subject. 4 

II. ANALYSIS 

We start with the coupled-line equations that describe the z evolu- 
tion of tA v (co, z) — the complex amplitude of the voltage in a mode v at 
position z resulting from a single Fourier component of the optical 
excitation at frequency co at position 2 = 0. (We can use Fourier 
components since the guide is linear in voltage.) 

dA ^' z) = £ A,( u , z)K v J(z) exp [i(i8, - ft,)*], (2) 

OZ n 

where <s4 v (u, z) (the complex amplitude in mode v) = A v (u, z) 
X exp [—*j9, (&>)«], 

0, = propagation constant for mode v (which is a function of 
co), and 
Knvf(z) = coupling coefficient between modes n and v (K*^ = — K VII ). 

If the input power excitation is a function of time, then the average 
power response at position z is given by the average of the square of the 
complex envelope, a„(z, t), at position z, i.e., 

(s v (z, 0) = (1 0,(2, t) | 2 ) = f(^v(o> + a, z)<A*{u, z)) exp [i(« + <r)Q 

X exp (—uat)dwd(w + a), (3) 

where the complex envelope a v (t, z) is the Fourier transform of <A X (a-, z) . 
We shall next derive eq. (1) by obtaining a differential equation for 
(d/dz)(s v (t,z)). 

2.1 Perturbation theory approach 

Following techniques used previously by Marcuse, 2 we first write 

f z (A v (u + a, Z)A*(<*, *)> = ^ [ f z A,(o> + a, *)] AtCa,, Z)) 

+ (A v (< a + * f z)[£ e A* v (u,z)]y- (4) 
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Using (2) in (4) we obtain 

— (A v (u + <r,z)A*(<a,z)) 

= <E A„(« + a, z)K v J(z) 
*» 

•exp {z'[/3„(" + a) - 0»(o> + a)-\z\A* v {o>, z)) 

+ <L ^->(" + *, *Wl(w, z)Kj(z) exp { -*rj3,(u) - /3„(o,)>}). (5) 

Next we recognize that the perturbation solution to (2) is given by 

A v (z) = A v (z') + f tf„ M A M ( 2 ') f f(x) exp [-H(/3, - M )x]d:c. (6) 

1 J z' 

We shall now substitute (6) into (5) with the following approxima- 
tions: We assume f(x) is independent of A v {z')Al(z') for x > z'. We 
keep only terms that are second order in K TI1 , because first-order terms 
would include the factor f(x) to first order and, since f(x) has been 
assumed independent of A v (z')A»(z') for x > z', the expectation of 
these terms would vanish [/(.r) has zero mean]. 

We obtain 

-p z (A v (o> + <r,z)A:(co,z)) 

= LI (A t (a> + a, z')A* v (u>, z'))K, s K v » [' </(*)/(*)> 

H 6 J V 

•exp iTj3 M (« + a) - Pt(o> + <t)](x — z)dx 
•exp {tTj3,(w + er) - 0a(« + <r)>} 

+ E E <A„(« + ^ fO^ifo z'))K*Xs f </(*)/(*)> 

/i 8 ■/*' 

•exp {— t|j8,(w) - 0a («)](* - z)}dx 

•exp {*rj8,(« + <r) - 0„(« + a) - 0,(») + &(«)>}, (7) 

plus two similar terms where w and co + a are interchanged and the 
conjugate is taken. 

In (7) we have again assumed independence between f(x) and 
-4 M (co, z') for x > z'. Now we assume that terms that are rapidly vary- 
ing in z can be neglected, i.e., we neglect terms proportional to 
exp{i[/3 M (w) - /3 (co)]2} unless n = v. We obtain for small z - z' [so 
that (A{u + a, z)A*(o>, z)) « (A (a; + er, g')A*(w, «'))]• 

4- {A.(u + a, «)A:(w, «)) = - E (A.(« + <r, 2)A:(co, 2))|/C m ,|»F 

+ E <ii M (« + a > 2M;(«,2)>|/i: M .|»F 

•exp {i[j8,(« + <r) - 0,(« + a) - /3,(a>) + /3 M («)>}, (8) 
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where F is the spectral height of the random process f(x), i.e., 

•exp {*'D3 M (« + o) ~ P.(« + *)](* ~ «)M(* ~ *) ( 9 ) 
is assumed independent of [?«(« + <r) — j8,(« + <r>], which implies 

(/(»)/(«)> «w(*-«). 

Next we assume that each mode has a well-defined group velocity 
within the frequency band of interest, that is, 

(S,( W + «r) -0,(")=^;- (10) 

Substituting (10) into (8) we obtain, using (2), 

a 

— (4,(u + a, z)A*(a, z)) 
az 

- ~ ««^,(« + <r, 2)c^(«, «)) exp (WO) 

02 

= - L |ff,,,| 2 F<«^,(w + <r, *)*/£(«, «)> exp (WO 

M 

+ E |X M .| S F<^„(« + v, z)<Al(u>, z)) exp (WO- (11) 
But 

a 

— (<c^ n (co + a, «)e^(w, s)> exp (WO) 
az 

= exp (WO — (c4 v (u + a, 2)^t(w, 2)) 

+ ^<^,(« + <r,a)^:(«,«)>]- (12) 

Substituting (12) into (11) and Fourier transforming as in (3), we 
obtain 

= - £ 1^.1 v<«-(*» *)> + E I ^.I'^a 2)), (13) 

which is the same as eq. (1). 

2.2 Kronecker product approach 

Rowe and Young 3 assume two modes and write (in our notation) 

dAl [ a ' z) = K l2 f(z) exp [»0fc - 0ffl A,(«, 2) 

33 (14) 

^'fo Z) = K n f(z) exp [t09, - /SO*] Ai(«, 2), 

02 
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where 

K\2 = K21 = iK. 

They then assume that f(z) is the derivative of an independent in- 
crements process satisfying 

</(z)/(z')> = F6(z - z') (Dirac delta). 

They divide the fiber into intervals (sections) of length A and assume 
that all the coupling takes place at the right end of each interval. They 
obtain the following difference equation : 

Ai(u, ZA) = cos (KcdAitu, (I - 1)A] + i sin (Kc^A^w, (Z - 1)A] 

•expp(/3i-/3 2 )ZA] 
A 2 («, ZA) = i sin (KcdAilu, (I - 1)A] exp [i(fc - 0OZA] 

-I- cos (Kci)A 2 [w, (I - 1)A], 

where 

ci = f(z)dz. 

J Interval I 

From (15) we obtain 

(Ai(u + <r, l&)Ai(o>, ZA)) 

= (cos 2 (Kci))(Ailu + <r, (Z - DAIKO, (I - 1)A]) 
+ (sin 2 (K Cl ))(A 2 lo> + <r, (I - l)A]^*[co, (I - 1)A]> 

•exp {t[/3i(o> + <r) - /3i(w) - /3 2 (o> + <r) + /3 2 (co)>}. (16) 

In (16) we used the fact that ci is independent of Ai[co, (Z — 1)A] 
and A 2 [u, (Z — 1)A] because f(z) has been assumed to be the deriva- 
tive of an independent increments process. We also used the fact that 
terms like (sin (Kci) cos (Kci)) equal zero because ci is a symmetrical 
random variable. 

In the limit as A — > (intervals { I } get small) we have 

(cos 2 (K C[ )) -+ 1 - FK 2 A (sin 2 Kci) -+ FK*A. 
We obtain from (16) 

— (Ai(oj + a, z)A* 1 (o), z)) 

- -Ftf 2 (Ai(o, + <x, z)Ai(<*, z)) + FK*(A 2 (o> + «r, 2)A 2 *(a>, 2)) 

•exp {i[/3i(<o + <r) - 0i(«) - /3 2 (a> + <r) + /3 2 (u)]z}, 

— <A 2 (a> + 0-, 2)A 2 (w, 2)) 

- -FK*(Atfa + «r)ilj(«)) + FiC 2 ^^ + <r)At(co)> 

•exp U'[/3 2 (o; + a) - /3 2 (co) - /3i(co + er) + /3i(co)>). 
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This is the same as eq. (8) with K n = Kn = iK and Ku = if 22 = 0. 
Thus, the coupled-power equations follow from (17) by Fourier trans- 
forming and making use of approximation (10). 
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